p 


AD-AO01  946 


6E0RGI A  INST  OF  TECH  ATLANTA  SCHOOL  OF  PHYSICS  F/S  7/4 

ELECTRONIC  EXCITATION  STRUCTURE*  CORRELATION*  AND  RELAXATION  AN— ETC(U) 
NOV  79  M  i  RIBARSKY  AFOSR-76-3106 

AFOSR-TR-80-0202  NL 


UNCLASSIFIED 


m  u. classified 


SECURITY  CLASSIFICATION  OF  THIS  PAPE  (Whit  DeleEnlered)^ 

~ - DOCUMENTATION  PAGE  beforeVomIle^Igkorm 

T73Ep3^  - ICSoVT  ACCESSION  NO.  *•  RECIPIENT'S  CAT  ALOC  NUMBER  * 


rs^. 


L  flTLE  (*fT Subtitle)  '  S.  TYPE  OF  REPORT  •  PERIOD  COVERED 

^Electronic  Excitation  Structure,  Correlation,  and  FINAL 

'Relaxation  and  Their  Effects  on  the  Optical  Oct.  1,  1976  to  Sept.  30 

Properties  and  Electron  Emission  of  Transition  /  6.  performing  org.  report  number 
Metal  Oxides^  - ...  -  / 


FINAL 

Oct.  1,  1976  to  Sept.  30, 


OKH7  _J _  -  - 

20  M.  W./Ribars^y^ 


^AFGSR- 7  6- 31^8  \ 


iUMBERfO 


9-  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS  _  *0.  PROGRAM  ELEMENT.  PROJECT.  TASK 

' \  _  AREA  A  WORK  UNIT  NUMBERS 

rzTJXX:  f  Iechnol°8y  [Y\)^L 

Atlanta,  Ga.  30332 _ — ~~ _ (JJgJ  _ 

It.  CONTROLLING  OFFICE  NAME  AND  AOORESS  12-  REPORT  DATE 

AIR  FORCE  OFFICE  OF  SCIENTIFIC  RESEARCH/NP  Nov.  30,  1979 

BOLLING  AFB,  BUILDING  410  '*•  number  of  pages 

WASHINGTON,  DC  20332  _ _ _ 

14.  MONITORING  AGENCY  NAME  a  ADDRESS?!/  dltlerent  from  Controlling  Olllee)  IS.  SECURITY  CLASS,  (ol  Ihle  rmpott) 


(mMi 


6,102^iiA5 

1*.  REPORT  DATE 

NOV.  30,  1979 

13.  NUMBEg^3F  PAGES 


(T)  ^iivdo  I  repti 
1  10bt  %-  £ 


unclassified 


tjyi 


15a.  DECLASSIFICATION/DOWNGRADING 
SCHEDULE 


I  16.  DISTRIBUTION  STATEMENT  (of  thla  Report) 


APPROVED  FOR  PUBLIC  RELEASE,  DISTRIBUTION  UNLIMITED 


HjUS 


$  ■  17.  DISTRIBUTION  STATEMENT  (ot  the  mbatract  entered  In  Block  20,  II  different  from  Report) 

H  (jJjZj?  Nov  11 


10.  SUPPLEMENTARY  NOTES 


1 19.  KEY  WORDS  ( Continue  on  ravaraa  aide  It  neceaaary  and  Identity  by  block  number) 


OPTICAL  PROPERTIES,  TRANSITION  METAL  OXIDES,  BAND  STRUCTURE,  CORRELATION 


20.  ABTTsPACT  (Cqnllnu*  on  raverae  aide  It  neceaaary  and  Idantlly  by  block  numbar) 

^Studies  of  the  electronic  excitation  structure  of  transition  metal  oxides 
and  the  effect  of  the  structure  on  the  optical  properties,  electron  correlation 
and  photoemission  spectra  of  these  materials  have  been  made.  Techniques  have 
been  developed  which  for  the  first  time  allow  precise  and  general  calculations 
of  the  excitation  structure  so  that  in  addition  to  the  above  properties  several 
other  properties  can  be  obtained.  These  techniques  have  been  applied  in 
particular  to  calculations  on  T10.  The  calculations  allow  comparison  with  . 


DD  ,^:M7,  1473 


EDITION  OF  I  NOV  *S  IS  OBSOLETE 


UNCLASSIFIED  jC  £  ^  h  ^ 


v  f  i  Afinr  «Tt 


sr  ViltE  r»B^*P*  'II'.  .n  *S.»,  M 


4  > 


mr- 


UNCLASSI FI  ED 


WCttfiTY  CLAHiriCATIQW  or  THIS  Mi  bi(«ra« 


-  2  - 


experiment  and  identification  of  the  structure  In  the  optical  spectra  of 
okidized  surfaces.  They  also  show  that  significant  correlation  and 
relaxation  effects  occur  in  the  band  structure. 


K 


i 


ELECTRONIC  EXCITATION  STRUCTURE,  CORRELATION 


AND  RELAXATION  AND  THEIR  EFFECTS  ON  THE 
OPTICAL  PROPERTIES  AND  ELECTRON  EMISSION 
OF  TRANSITION  METAL  OXIDES 


P/'H  (A>'SC) 

■d  r\nd  Is 

-^-12  17b). 


Final  Scientific  Report 


Electronic  Excitation  Structure,  Correlation,  and  Relaxation 
and  their  Effects  on  the  Optical  Properties  and  Electron 
Emission  of  Transition  Metal  Oxides 

By 

M.  W.  Ribarsky 

November  30,  1979 


Grand  No.  AF0SR-76-3108 


TABLE  OF  CONTENTS 


Page  No. 

Introduction .  1 

Objectives . . .  3 

Achievements 

A.  Development  of  Dielectric  Response  Theory  in  Local  Orbitals 

Framework.. . . . . . . . .  4 

B.  Devising  Methods  of  Calculation...... .  8 

C.  Investigation  of  Transition  Metal  Oxides .  14 

References .  27 

Communications .  28 

Personnel . 29 

Appendix 


■— i ww  am  .  Htiommm 


I.  INTRODUCTION 


A  system  responds  to  an  external  probe  by  undergoing  various  excitations 
In  order  to  distribute  the  energy  and  momentum  transferred  from  the  probe.  The 
nature  of  these  excitations,  their  characteristic  energies  and  distributions, 
establishes  the  properties  of  the  system.  Thus  the  electronic  excitation  struc¬ 
ture  of  a  system  establishes  its  optical  and  electrical  properties  and  the 
pattern  of  its  electronic  emission. 

An  understanding  of  the  electronic  excitation  structure  elucidates 
experimental  methods  which  probe  the  optical  properties,  the  electron  emission 
and  energy  loss  spectra  and  the  collective  excitation  modes  of  a  system.  Ad¬ 
vancing  in  an  interlocking  manner,  the  experimental  methods,  which  are  now 
providing  a  wealth  of  precise  data,  then  give  an  even  more  refined  interpretation 
of  the  theory  which  can  be  revised  and  tested  anew.  Moreover,  when  the  relation¬ 
ship  between  experiment  and  theory  is  understood,  experiment  can  be  used  to  show 
the  effect  of  changes  in  the  basic  electronic  properties  and  other  properties 
not  included  in  the  theory.  For  example,  the  probe  of  oxides  grown  an  titanium 
with  reflectivity  and  Auger  measurements  yields  information  on  the  basic  electronic 
structure  of  the  oxide  which  can  be  related  to  changes  which  also  occur  in  the 
electrical  conductivity  and  growth  rate  of  the  oxide. 

In  this  work  a  dielectric  response  theory  has  been  formulated  to  describe  the 
electronic  excitation  structure  of  a  solid  state  system.  The  properties  and 
experimental  methods  mentioned  above  can  be  investigated  within  the  framework  of 
this  theory.  In  addition  one  can  look  at  electronic  correlations  which  can  produce 
large  collective  modes  such  as  plasmons  or  which  can  have  profound  effects  on  the 
band  structure  of  the  system. 

We  have  used  the  local  orbitals  technique  in  constructing  the  dielectric 
response.  This  framework  allows  the  wedding  of  ground  state  calculations  of  energy 
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bands  and  wavefunction  with  excitation  calculation  in  the  same  representation, 
and  it  is  especially  well  suited  for  materials  with  narrow  d-bands.  With  the 
local  orbitals  technique,  we  can  calculate  symmetry-breaking,  localized  exci¬ 
tations  more  easily.  The  dielectric  response  model  then  allows  the  inclusion 
of  both  the  symmetric-state  ground  and  the  symmetry-broken  excited  state  and 
permits  the  calculation  of  the  excitation  energy  and  transition  rate  in  one 
step. 

We  have  concentrated  on  the  properties  of  transition  metal  oxides  in  this 
work.  The  partially  filled  d-bands  are  the  distinguishing  feature  of  transition 
metal  chemistry.  These  bands  are  relatively  narrow,  but  they  are  overlapped  by 
the  next  higher  s-band  causing  a  great  deal  of  s-d  hybridization.  On  the  other 
hand  in  some  transition  metal  compounds,  such  as  many  of  the  oxides,  the  bonds 
are  mostly  ionic  and  the  highest  metallic  s-band  is  unoccupied.  In  this  case 
the  highest  s-band  moves  away  from  the  partially  filled  d-band.  Then  the  s-d 
hybridization  is  much  less  and  the  d-band  is  quite  narrow.  Quantitative  know¬ 
ledge  of  the  degree  of  hybridization  of  these  bands  and  the  related  structure  of 
the  bonds  in  the  solid  is  very  important  for  understanding  many  mechanisms  such 
as  electrical  conductivity  and  catalysis.  Since  many  of  the  oxides  of  the  trans¬ 
ition  metals,  especially  titanium,  exist  in  several  structures  with  widely  dif¬ 
ferent  electrical  properties,  this  knowledge  also  enables  us  to  identify  the 
oxides  occuring  on  transition  metals  by  analyzing,  for  example,  optical  absorption 
spectra.  Evidently,  the  electronic  structure  of  these  oxides  is  related  to  the 
rate  of  oxidation  of  the  transition  metal  and  to  its  degree  of  corrosion. 
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II.  OBJECTIVES 


The  purpose  of  this  research  was  to  formulate  and  implement  a  dielectric 
response  method  for  calculating  excitation  structure  and  to  apply  the  method  to 
calculations  of  the  electronic  properties  of  transition  metal  oxides.  Important 
initial  objectives  were  to  set  up  a  dielectric  response  theory  in  the  framework 
of  the  local  orbitals  technique  and  to  devise  methods  for  calculating  the 
excitation  structure.  This  latter  objective  was  not  easily  realized  since  it 
required  very  efficient  techniques  to  handle  large  matrices  and  to  calculate  and 
sort  large  numbers  of  multicenter  two  electron  integrals.  The  next  objective 
was  to  focus  on  transition  metal  oxides,  especially  TiO.  The  properties  of  these 
materials  were  analyzed  with  emphasis  on  the  optical  properties.  Experimental 
results  were  viewed  in  the  light  of  this  analysis,  and  several  properties  of  the 
oxide  system  under  different  conditions  of  oxidation  were  elucidated.  A  final 
objective  was  the  study  of  several  materials  to  get  an  idea  of  the  range  of 
excitation  structures  and  the  wide  variation  of  electronic  properties. 
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III.  ACHIEVEMENTS 


A.  Development  of  Dielectric  Response  Theory  In  Local  Orbitals  Framework. 

The  central  point  of  any  local  orbitals  expansion  is  that  the  charge  density 
may  be  written  as 


p(r,r') 


l 

ij 

a,b 


Tai,bj  Vr’> 


(1) 


where 

<f>ai(r)  =  ^(r-Ra), 


and  <|>  is  localized  in  space.  The  matrix  T  contains  the  coefficients  for  the 
Bloch  states  in  terms  of  linear  combinations  of  the  <j>’s.  A  particular  choice  of 
<J>  occurs  when  one  rotates  in  the  subspace  of  occupied  Bloch  states.  This 
rotation  leaves  p  in  Eq.  (1)  unchanged,  but  allows  one  to  solve  for  it  in  a  new 
set  of  equations 

(F(p)  +  pAp)|*ai>  -  Eai|*ai>  (2) 


where  F  is  the  Fock  operator  and  A  is  an  arbitrary  Hermitiar.  operator  chosen 

such  that  <(>  .  is  well- localized.  This  is  the  "local  orbitals  methods"  of  Kunz1, 
ai 

the  results  of  which  we  have  used  extensively  in  this  work. 

One  may  use  the  field-theoretic  methods  to  contruct  the  particle  hole  (PH) 
polarization  propagator  Mo>),  the  poles  of  which  are  the  excitation  energies. 

Under  certain  assumptions, 2  the  equation  for  IT  (to)  reduces  to  a  Dyson-like  equation 

n(w)_1  C(w)  =  A(u>)  C(u>)  (3) 


with 


n  (m)-1  «  n°(u!)-1  -  K(u>) 


(4) 
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The  kernel  K(<d)  contains  two  electron  matrix  elements  describing  the  particle- 

hole  interactions.  n_°(to)  is  the  zero' th  order  PH  propagator  (without  particle- 

hole  interactions)  and  thus  is  constructed  wholly  with  ground  state  energy  bands 

and  wave  functions.  One  solves  Eq.  (3)  for  the  values  id  =  to  for  which  A(id)=0. 

n 

These  are  the  excitation  energies  of  the  crystal.  Then  n (id)  is 

c^(w)c  ,  (w)+  c"(<d)c>)  + 

n  .(<*>)-  X  — — 7-frj - 5 -  ,c, 

ss  L  id— u )  (w)+in  id-id  (w)-iri  (5) 

n  n  J 

n 

+  2 
where  n  =  0  .  In  the  random  phase  approximation  (RPA) ,  K(id)  does  not  depend 

on  id  so  that  and  are  independent  of  id.  Then  the  solution  of  Eq.  (3)  is 

much  simplified  as  we  will  show  below. 

The  difficulty  in  solving  Eq.  (6)  for  crystalline  systems  is  that  II_((d)-1 
is  not  block  diagonal  in  momentum  as  is  the  Fock  operator,  for  example.  It 
depends  on  the  initial  momentum  1c  and  on  the  momentum  transfer  between  ground  and 
excited  state  q.  Thus  11^  has  off-diagonal  elements  in  momentum  which  give  rise  to 
"local  field  effects".  These  effects  are  related  to  the  formation  of  a  "localized 
hole"  upon  excitation,  for  if  the  hole  were  delocalized,  the  excited  state  would 
have  the  full  symmetry  of  the  ground  state  and  n.(u>)  would  be  diagonal.  Thus  when 
local  field  effects  are  present,  one  must  diagonalize  a  huge  matrix  in  momentum 
space.  On  the  other  hand,  if  instead  of  momentum  a  sufficiently  localized 
representation  such  as  Wannier  functions  were  used,  one  could  diagonalize  n_(oj)— 1 
in  a  manageable  subspace. 

We  can  avoid  the  difficulties  of  the  Bloch  representation  by  transforming  to 
a  local  orbitals  basis.  The  transformation  requires  a  matrix  Tg(k1k2)  which 
contains  all  the  coefficients  for  a  local  orbitals  expansion  of  the  Bloch  functions 
where  S  is  a  condensed  notation  for  the  local  orbitals  indices,  and  kj/k2  are 
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are  condensed  notations  for  the  initial/final  Bloch  states.  Then 

n°  =  l  t  *(k1,"k2+<i)  x°Cki,k2+q  ;w)-1  t  , (ki,tt2+q)  (6) 

ss  k!k2 

Xq-1  is  just  the  inverse  of  the  polarizability  for  the  non-interacting  electron- 
hole  pair3,  and  it  is  given  by 


-  vi-W  ‘“'V'5  +  V  <7> 

-  {" +  %+;  -  v 

where  E*  and  F>  are  the  one-particle  energy  and  occupation  number  for  Bloch  state 
1c.  The  inclusion  of  K  in  Eq.  (4)  introduces  particle-hole  interaction  effects. 

In  the  RPA,  K  is 


K  .  v„.i  ui4t„».,  ?  ~  „i,.i  k'n 


ss 


ai  bj,  b'j'a'i'  ai  a'i  ,b  j'bj 

fr 


aibj.b  j  a  i  ai  b  j  .  . 

fr-r' 


bj(r)<?a'i'(r’)V- 


(8) 


The  set  of  equations  (3-8)  constitute  a  system  to  find  the  energies  and 
transition  rates  for  excitations  from  the  self-consistent  ground  state  of  the 
crystalline  solid.  A  much  smaller  subspace  of  local  orbitals  is  needed  for  these 
equations  than  would  be  necessary  using  Bloch  functions.  In  addition  the  calculation 
is  based  on  the  fully  self-consistent  ground  state  band  structure  of  the  system- 
this  is  the  first  time  such  a  calculation  has  been  done.  Hanke  and  Sham4  have 
used  a  similar  set  of  equations  in  model  calculations  on  diamond.  However  they 
directly  invert  n_((u)_1  in  Eq.  (4)  for  each  ui;  such  an  approach  would  be 
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prohibitively  time-consuming  for  a  full  calculation.  As  we  will  show  below,  it 
is  possible  to  diagonalize  Eq.  (3)  and  get  a  complete  spectrum  in  one  step. 

To  proceed  we  may  decompose  Eq.  (6)  using  Eq.  (7)  thusly 


n°s,(q,w)-1  =  u  0ss,(q)  Wss,(q) 


(9) 


where 


0  '<q)-I  Ta *(kx,k2+q)  (F.  T  ,(t1(t;2+g)  (10) 


ss  ■  -  f  s 
kik2 


WSS’(9)  =  l  Tg  (ki,t2+q) 
kjk2 


kj  k2+q  s 


(F>  +  F>  -  2P>  F*-  >) 

k2+q  k1%-Kj; 


(11) 


x  (E>  +  -  E>  ) 

k2+q  k2 


Tsi (kltk2+q). 


0(q)  in  Eq.  (10)  behaves  like  an  overlap  matrix  in  that  it  is  positive  definite. 
Then  we  can  write 


0(q)  =  L(q)L+(q) 


(12) 


so  that  for  A  =  0  and  after  transformation,  Eq.  (3)  becomes 


Wn(q)Cn  (q)  =  [W' (q)  -  K' (q) 


(q) 


(13) 


with 

W' (q)  =  L~ 1 (q)W(q ) (L+(q)" 1 

and  so  on.  The  nature  of  the  transformation  ensures  that  the  matrices  in  Eq.  (13) 
remain  Hermitian.  We  can  therefore  treat  Eq.  (13)  as  an  eigenvalue  problem  and 
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diagonalize  to  find  the  spectrum  for  each  q.  For  the  optical  spectrum,  we 
may  treat  the  limit  q+0  and  diagonalize  Eq.  (13)  only  once.  However,  for 
problem  such  as  correlation,  a  smapling  for  q  throughout  the  Brillouin  Zone 
must  be  made. 

B.  Devising  Methods  of  Calculation. 

It  turns  out  that  the  solution  of  Eq.  (13)  is  no  easy  matter  because  of  the 
size  and  complexity  of  the  matrices  involved  and  because  of  the  need  to  calculate 
multicenter  integrals.  We  discuss  some  of  the  methods  used  on  this  problem  in 
more  detail  in  this  section  because  they  represent  achievements  of  wider  applica¬ 
bility  than  just  to  the  present  problem. 

As  presented  in  Eq.  (8),  the  kernel  K  depends  on  individual  multicenter  two 
electron  integrals.  Two  problems  exist  with  the  calculation  of  K.  The  first  is 
that  the  multicenter  integrals  can  be  difficult  and  time-consuming  to  calculate. 

In  fact,  if  one  uses  Slater  type  orbitals  (STO's)  the  integral  expression  does 
not  have  a  closed  form  and  one  must  do  an  expansion  over  angular  momentum  states, 
the  sum  becoming  more  complicated  for  three  -  and  four  -  center  integrals.  The 
second  problem  is  that  the  integrals  do  not  sum  nicely  into  a  quantity,  like  the 
Fock  matrix  in  Hartree  Fock  calculations,  which  has  the  initial  symmetry  of  the 
system.  This  happens  because  the  excited  system  no  longer  depends  on  just  the 
one-electron  density  matrix  and  completely  breaks  symmetry  in  general.  Thus  the 
sorting  of  the  integral  labels  is  a  major  accounting  problem.  We  have  found  quite 
effective  methods  for  solving  both  these  problems  in  the  Unrestricted  Hartree 
Fock  (UHF)  molecular  cluster  codes  developed  by  Robert  Euwema  and  Frank  Tobin 
under  an  Air  Force  grant.  Frank  Tobin  has  consulted  with  us  and  helped  us  to 
understand  the  methods  used  in  the  UHF  code.  He  has  also  helped  us  to  maintain 
the  code  in  running  order,  and  we  are  using  it  for  several  other  calculations. 
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Initially,  since  our  original  local  orbitals  were  in  STO's,  we  calculated 
our  multicenter  integrals  in  terms  of  STO's.  We  devised  a  set  of  codes  which 
would  calculate  integrals  for  states  of  any  angular  momentum  on  any  set  of  sites. 
These  codes  will  be  useful  for  limited  calculations  in  the  future  (for  example, 
calculations  involving  high  angular  momentum  states),  but  they  proved  too  slow 
for  the  excitation  calculations.  However,  the  UHF  code  has  the  most  efficient 
multicenter  integrals  package  currently  available  for  calculations  involving  s,p, 
or  d  symmetries,  and  we  were  able  to  modify  the  package  for  our  purposes.  This 
code  calculates  integrals  in  terms  of  Gaussian  type  orbitals  (GTO's),  and  the 
integrals  can  be  calculated  without  an  expansion  in  terms  of  angular  momentum 
states.  The  extra  effort  necessary  to  expand  the  individual  STO's  in  terms  of 
GTO's  was  worth  paying  in  return  for  the  much  quicker  integral  generation. 

The  sorting  of  the  integrals  labels  was  complicated  by  the  fact  that  some  of 
the  lattice  vector  indices  must  be  summed  over,  for  example,  in  the  optical 
absorption  calculation  before  the  final  sort  can  be  done.  Luckily,  the  UHF  package 
has  a  remarkable  general  code  based  on  the  Yoshimine  procedure  which  allows  the 
sorting  of  labels  initially  in  any  order.  We  were  able  to  adapt  this  code  to  our 
purposes. 

After  the  matrices  are  finally  constructed,  we  face  the  problem  of  diagonal¬ 
izing  them.  Actually,  two  diagonalizations  must  be  done  on  large  matrices  of 
somewhat  different  character.  The  matrices  are  large  because  two  sets  of  indices 
are  needed,  one  set  for  the  initial  electron  and  one  for  the  final  state,  and 
because  the  lattice  vector  must  also  be  Indexed.  In  general,  all  the  indices 
appearing  on  the  RHS  of  Eq.  (8)  explicitly  appear  in  the  matrix.  For  our 
calculation  on  TiO,  where  it  was  only  necessary  to  use  the  valence  p  and  d  local 
orbitals  and  up  to  the  next  nearest  neighbor  lattice  vectors,  the  dimensions  of 
the  matrices  were  still  622  x  622.  For  a  complex,  Hermitian  matrix  of  this  size. 
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standard  diagonalization  techniques  would  take  about  three  hours  of  computer  time 
to  find  all  the  eigenvalues  and  eigenvectors.  We  devised  a  scheme  involving 
diagonalization  of  sub-blocks  within  the  matrix  and  then  correcting  them  pertur- 
batively  which  markedly  reduced  the  computation  time. 

The  matrix  C)  in  Eq.  (10)  is  dense  and  does  not  easily  separate  into  blocks. 
However  it  is  diagonally  dominant  and  blocking  tendencies  are  apparent.  To 
illustrate  the  procedures  used,  in  Fig.  1  we  show  density  maps  of  a  150  x  150  sub¬ 
block  of  our  large  matrix.  Each  symbol  in  the  maps  represents  an  average  magnitude 
of  the  elements  in  the  small  block  at  that  position,  the  larger  magnitude  elements 
being  represented  by  darker,  built-up  symbols.  The  blank  areas  in  the  maps 
contain  non-zero  elements  which  are  below  a  given  tolerance  in  magnitude.  Thus 
the  map  is  a  display  of  the  distribution  of  magnitudes  of  matrix  elements.  The 
first  step  is  to  take  the  original  matrix  in  Fig.  la  and  interchange  rows  and 
columms  so  that  the  diagonal  elements  are  ordered  in  magnitude.  As  shown  in 
Fig  lb,  this  step  greatly  condenses  the  density  pattern.  We  then  choose  blocks 
along  the  diagonal  (in  this  case  four  are  chosen)  so  that  the  natural  blocking 
pattern  displayed  by  the  map  is  followed,  and  we  find  the  eigenvalues  and 
eigenvectors  of  each  block.  We  perform  the  unitary  transformation  using  the 
block  diagonal  matrix  composed  of  the  eigenvectors  with  the  result  in  Fig.  lc. 

Now  the  diagonal  blocks  just  have  eigenvalues  in  descending  order  on  the  diagonal 
with  zeros  on  either  side,  and  the  off-diagonal  blocks  also  have  a  particular 
order.  The  matrix  is  now  ready  for  further  rearrangement.  One  possiblility  is 
shown  in  Fig  Id  where  we  have  simply  sorted  the  diagonal  elements  again  as  we 
did  in  Fig.  lb.  The  solutions  are  finally  corrected  perturbatively  as  outlined 
below.  Before  the  perturbation  corrections  is  attempted,  however,  eigenvalues 
which  are  very  close  or  which  are  connected  by  large  off-diagonal  matrix  elements 
must  be  blocked  together  since  they  will  make  the  perturbation  expansion  unstable. 
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For  optimally  blocked  matrices  two  diagonalizations  of  sub-blocks  followed  by 

perturbative  corrections  give  good  results.  For  diagonalizing  any  large,  non-sparse 

matrix,  a  topic  which  has  not  received  much  attention,  these  procedures  are  useful. 

The  matrix  W' (q)-K' (q)  in  Eq.  (13)  is  also  dense,  but  it  blocks  more  easily  owing 

to  the  contributions  of  the  two  electron  terms  in  K' (q).  We  can  diagonalize  each  of 

these  blocks  separately,  and  then,  using  a  variation  on  degenerate  perturbation 

theory,  we  can  include  the  effects  of  terms  not  in  the  blocks.  Let  us  assume 

we  want  to  block  diagonalize  the  matrix  n_(uj)-1  which  is  arranged  as  in  Fig.  2. 

We  have  assumed  the  total  matrix  is  of  size  Nm  x  Nm  and  have  arranged  the 

elements  along  the  diagonal  into  Nmxm  blocks,  n-1 . IT-1,  although  there  is  no 

— 1 1  — nn 

need  to  restrict  ourselves  to  blocks  of  equal  size.  Then  the  eigenvalue  problem 
becomes 


A-n  )c  -  l  n“1(A-n"1)  xn 
~11  l  ~i„  “wm  -ui  i 


-l  -l  -l  -l  -l  -l  -i_ 

+  l  n  (A-n  )  n  (A-n  )  n  _c,+.... 

In  “UP  ~PVV  ”VV  — V 1  l  * 

P^v^l 

where  A  is  a  true  eigenvalue  of  the  total  matrix  and  is  a  projection  of  the 

-1 

total  eigenvector  on  the  subspace.  The  elements  of  could  be  arranged  in  order 

-1  -1  -l 

of  importance  with  respect  to  so  that,  for  example,  H  is  larger  than  n  and 

11  12  “13 

so  on.  In  practice,  the  series  in  Eq.  (22)  is  rapidly  convergent;  and  if  one  uses 

some  of  the  optimization  procedures  outlined  above,  the  first  term  on  the  RHS  of 

Eq.  (14)  is  all  that  is  necessary  for  good  results.  Note  that  if  we  transform 
-1  -1  “1 

lL(w)  as  we  did  in  Fig.  lb,  (A-n_  )  is  just  the  inverse  of  a  diagonal  matrix. 


Figure  2.  Block  Form  for  PH  Propagator  Matrix. 
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C.  Investigation  of  Transition  Metal  Oxides. 

1.  Optical  Properties  -  These  investigations  center  especially  on  the 
properties  of  TiO.  Our  calculations  and  analyses  of  experiments,  indicate  that 
the  optical  reflectivities  of  titanium  surfaces  undergo  significant  changes  as 
oxides  are  formed.  These  calculations  identify  previously  unidentified  peaks 
in  the  oxide  spectra  at  about  8-10  eV  and  26  eV  as  transitions  from  oxygen  2p 
and  2s  bands,  respectively.  Our  calculations  also  identify  structures  in  the 
5-6  eV  region  as  involving  d-bands.  The  d-band  structure  changes  quite  a  lot 
under  varying  conditions  of  oxidation  and  can  be  used  as  an  indicator  of  the 
oxide  structure. 

Our  analysis  indicates  that  the  oxide  produced  initially  by  oxygen  bombardment 
of  the  clean  polycrystalline  surface  is  TiO-like5  with  further  banbardment  producing  a 
trend  toward  greater  oxidation  or  perhaps  a  mixed  oxide  state  but  that  the  air  exposed 
initial  surface  seems  to  be  Ti02«  Oxygen  exposure  of  the  clean  surface  at  room 
temperature  p  roduces  a  TiO  structure  which  evolves  at  a  much  slower  rate  than  for 
the  single  crystal.6  The  changes  in  the  reflectivity  data  which  accompany  these 
different  conditions  of  oxygen  exposure  are  consistent  with  the  assumed  oxide 
structures.  For  example,  in  Fig.  3  we  show  the  changes  in  the  reflectivity  which 
accompany  continued  oxygen  bombardment.  The  structure  around  5  eV,  identified  by 
our  calculations  as  involving  d-bands,  reduces  dramatically  which  is  consistent 
with  a  trend  toward  Ti02,  in  which  case  the  valence  d-bands  would  be  completely 
depleted.  Undoubtedly  the  condition  of  oxygen  bombardment  in  these  measurements 
prevents  the  surface  from  stabilizing.  Certainly  the  differences  among  the  oxides 
would  be  quite  important  in  the  mechanisms  of  oxide  growth  under  various  conditions, 
and  these  changes  may  well  be  useful  in  understanding  the  corrosion  and  possivation 
of  metallic  surfaces. 
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Let  us  look  in  more  detail  at  the  analysis  of  the  TiO  optical  properties. 

The  major  components  in  the  optical  absorption  spectra  of  all  the  oxides  are 
several  structures  due  to  a  d-band  excitations  at  lower  energies,  a  broader 
structure  due  mostly  to  0-2p  band  excitations  (perhaps  with  some  admixed  d 
character)  at  intermediate  energies,  and  a  structure  due  to  0-2s  excitations 
at  higher  energies.  All  these  interband  structures  are  superimposed  on  a 
temperature-dependent  intraband  background  which  has  a  Drude-like  frequency 
dependence  of  uT3  at  lower  frequencies.  The  positions  and  strengths  of  the 
interband  peaks  vary  with  varying  oxide  structure  so  that  the  optical  data 
provides  a  means  for  identifying  the  oxide.  As  an  example  of  the  optical 
behavior,  we  present  the  reflectance  spectra5  of  various  Ti  surfaces  in  Fig.  4. 

The  lower  two  curves  are  for  oxide  structures  formed  under  different  conditions. 
These  curves  show  peak  structures  from  about  5  to  6  eV  to  lower  energies  for 
which  d-band  transitions  are  assigned  (actually,  recent  low  temperature  measure¬ 
ments  show  the  d-band  structure  extends  down  to  .2  eV.  7).  Transitions  from  the 
0-2p  band  region  are  assigned  to  the  broad  structure  around  10  eV,  which  makes 
sense  since  this  structure  disappears  for  the  clean  surface.  The  assignment  for 
the  structure  seen  above  20  eV  which  peaks  somewhere  around  25  eV  is  0-2s,  although 
transitions  from  the  valence  bands  to  higher  conduction  bands  are  also  present. 

To  investigate  this  optical  behavior,  we  need  to  look  at  the  imaginary  part 
of  the  macroscopic  dielectric  function 


e(w)  =  lim  e(q,q;u>) 
-+ 

q*0 


Im(e(w))  may  be  written8 

Im(e(w))%-^  l  F®  Imnss,  (u>)  f“,’ 


(O'-  “  ,  S 

S,  S 


Fa  =  <;  Jv  | <t>  ,> 

s  ai1  a1  oj 


(15) 


(16) 


Reflectance  of  clean  and  oxidized  Ti  surfaces 


where  the  reducible  polarization  matrix  n_  has  the  structure  given  in  Eq.  (5). 

(X 

One  may  neglect  the  effect  of  the  matrix  elements  F^  in  Eq.  (16)  and  the 
particle-hole  interactions  contained  in  I[.  Under  these  approximations,  the 
macroscopic  dielectric  function  is  proportional  just  to  the  joint  density  of 
states. 

To  evaluate  the  importance  of  different  effects  on  the  optical  behavior,  we 
present  in  Fig.  5  the  joint  density  of  states  calculated  by  us  using  uncorrelated 
HF  energy  bands. ^  For  comparison  we  present  in  Fig.  6  the  experimental  joint 
density  of  states  obtained  from  the  Kramers-Kronig  analyzed  optical  constants  via 
the  relation 


Nt(oj)  =  id2  e,  (u) 

J  exp  2  exp  (17) 

where 

e(w)  =  ti(u))  -  iC2(to) 

We  see  that  there  is  general  agreement  among  the  positions  of  the  peak 
structures  in  the  two  spectra  as  indicated  in  Table  1.  In  order  to  identify  the 
symmetry  of  the  transitions  involved  we  may  consider  the  function 

iV- (u)  =  l  |cfe  vi(l0|2  5  (w-Eb(k) )  (18) 

b,1c 

where  c,  . (k)  is  the  Bloch  coefficient  for  local  orbital  i  in  band  b.  Thus  n  (»>) 
b*  l  i 

is  a  measure  of  the  density  of  states  for  a  given  local  orbital.  It  turns  out  that 
orbitals  of  the  same  symmetry  have  local  densities  of  states  which  are  quite  similar 
in  the  valence  and  conduction  band  regions  so  that  we  may  plot  representative  0-2p 
and  Ti-3d  orbitals  in  Fig.  7.  Actually,  the  orbital  densities  of  states  plotted  are 
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Fig.  6  Experimental  joint  density  of  states. 


Table  I.  Comparison  of  peak  maxima  from  the  calculated  and  experimentally 


derived  joint  densities  of  states  and  Lynch's  £3  data.  All  energies 
are  in  eV. 


a 

Calculated 

b 

Experiment 

c 

Lynch  et.  al 

3.28 

3.32 

3.2 

5.44 

5.20 

5.6 

6.34 

6. 45 (shoulder) 

9.40  (central  maximum) 

8.78 

8.2 

31. 

24.5 

23.8 

a.  See  Figure  5 


b.  See  Figure  6 


N=  ENERGY! RY I 

Local  densities  of  states  in  arbitrary  units,  a)  Contribution  of  the 
the  Ti-3d  local  orbital  in  the  valence  and  conduction  band  regions, 
b)  Contribution  of  the  0-2p  local  orbital. 
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quite  a  bit  larger  than  any  of  the  others  in  the  valence  band  region  and  are 
also  about  the  same  size  as  each  other,  the  differences  in  Fig.  7  being 
attributable  to  different  scalings.  Thus  the  data  in  Fig.  7  supports  assignments 
of  p  and  d  band  regions  in  the  optical  spectrum,  and  results  not  shown  here  support 
our  assignment  of  the  peak  at  about  25  eV  to  0-2s  transitions.  Note  that  n^(iu) 
in  Fig.  7  is  most  prominent  above  the  Fermi  energy  at  -.888  eV,  but  that  there  is 
some  mixing  between  the  p  and  d  states  in  the  region  just  below  the  Fermi  level. 
Whether  this  mixing  actually  occurs  in  the  real  system  will  be  discussed  below. 

The  first  effect  we  may  consider  to  improve  our  calculation  is  to  include 
the  matrix  elements  F^  in  Eq.  (16).  These  matrix  elements  have  the  effect  of 
optical  selection  rules  since  they  only  connect  states  with  angular  momentum 
difference  A£=+l.  This  effect  can  be  quite  significant  for  transitions  between 
bands  each  of  a  particular  symmetry,  but  the  effect  would  be  negligible  for 
transitions  between  well-mixed  bands.  To  see  what  is  the  case  for  TiO,  we  have 
recalculated  the  optical  absorption  with  the  matrix  elements  F™  included  and  have 
presented  the  results  in  Fig.  8.  The  shape  is  nearly  the  same  as  in  Fig.  5, 
indicating  the  symmetry  admixing  depicted  in  Fig.  7,  although  there  is  a  slight 
reduction  of  the  structure  at  lower  energy  because  the  selection  rule  forbids  d- 
to-d  transitions,  which  are  stronger  in  this  region.  Relaxation  and  correlation 
effects,  discussed  below,  will  separate  the  p  and  d-bands  so  that  the  lower 
energy  structure  will  be  further  diminished. 

2.  Relaxation  and  Correlation  Effects.  To  go  beyond  the  HF  approximation  for 
the  ground  state  properties  of  the  system,  it  is  necessary  to  consider  self  energy 
corrections.  These  corrections  introduce  relaxation  and  correlation  effects  to 
the  single  electron  ionizations  described  by  the  energy  band  structure.  The 
eigenvalue  equation  may  be  written 


J!o  +  E  (w) 


A(<d) 


E(w)  A (ui) 


(19) 
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Calculated  joint  density  of  states  with  matrix  elements  F  included 
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where  HQ  is  the  HF  Hamiltonian,  £  is  the  self-energy,  and  E  is  our  correlated 
band  energy.  Neglecting  possible  spin  polarization  effects,  we  may  write  the 
most  important  term  in  E  as  a  dynamically  screened  exchange10 

llk,k'(a>)  =  TiT  ^ki  Vi£,mn  Vm'n',£j 


x  /  IT  ,  ,  (q,u>* )  G°  ,  ,  .-+  | .  i6w'  , 

mn,m  n  kS.,  k  H  (q,w+u>  )  e  dw 


(20) 


where  k=kn,  k' =kn'  (The  reduced  momenta  are  the  same).  To  avoid  even  more 
complicated  notation,  repeated  indices  in  the  equations  for  £  indicate  that  the 
terms  are  to  be  summed  over  that  index.  In  Eq.  (20)  =  a£.  e  J (no  sum) 

where  a>^  is  the  Bloch  coefficient  for  state  k.  G°  is  the  one-particle  Green's 
function  constructed  with  HF  energy  bands,  and  J1  is  given  in  Eq.  (5)  (Except  that 
now  it  has  an  explicit  q  dependence).  One  can  easily  carry  out  the  w'  integration; 
the  result  contains  contributions  from  the  poles  of  n  and  from  the  poles  of  G°.10 

A  difficulty  with  Eq.  (20)  is  that,  unlike  the  optical  absorption  calculation 


where  only  q>o  is  needed,  one  must  find  JHq,w)  for  a  sampling  of  q  throughout  the 

Brillouin  zone.  In  our  model  this  requires  a  repeated  series  of  diagonalizations 

of  large  matrices.  An  alternative  which  retains  some  of  the  major  effects  we  seek 

sc 

but  is  computationally  much  simpler  is  to  use  the  second  order  expression  for  Z 
We  would  then  replace  FI  by  in  Eq.  (6).  One  can  show  that  such  an  approximation 
neglects  effects  due  to  multiple  particle-hole  excitations  or  collective  modes  such 
as  plasmons,  but  studies  on  highly  inhomogeneous  atomic  systems10  for  which  the 
electron-electron  interactions  should  be  strong  indicate  the  correct  relative  sizes 
for  relaxation  and  correlation  contributions  are  given  by  the  second-order  calcu¬ 
lation.  For  solid  state  systems,  where  collective  modes  are  more  significant,  the 
second-order  calculations  would  not  describe  regions  of  the  photoemission  spectrum, 
for  example,  where  plasmons  dominate.  However,  one  can  immediately  write  down  the 
expression  for  n_°  without  doing  any  diagonalizations  or  inversions,  and  the  two- 
electron  integral  sums  are  similar  to  the  ones  we  did  for  the  optical  absorption  case. 
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The  first  step  in  the  calculation  is  to  transform  Eq.  (19)  to  the  basis  of 


HF  Bloch  states.  Then  we  have 


E  I  +  I(w)  A(w)  *  E(u>)  A(w),  (21) 

o  — 

where  J.  is  the  identity  matrix.  The  effect  of  Eq.  (21)  is  to  admix  the  HF  band 
states.  But  u  appears  in  the  equation  as  a  parameter,  and  a  complete  solution 
for  the  band  structure  could  be  written  schematically  as 

E  =  E  +  X(E). 
o 

Therefore  we  want  solutions  of  Eq.  (21)  for  values  of  w  near  the  band  energies  of 
interest.  In  practice  we  solve  Eq.  (21)  for  a  sampling  of  frequencies  in  the 
energy  range  (e.g.,  the  valence  band  region)  we  are  studying. 

We  have  applied  the  above  method  to  calculating  relaxation  and  correlation 
corrections11  to  the  HF  energy  bands  of  Jennison  and  Kunz.9  Originally  these 
bands  were  thought  to  be  in  agreement  with  XPS  measurements12  of  the  0-2p,  Ti-3d 
valence  band  density  of  states  for  TiO,  but  subsequent  measurements13  showed  the 
first  experiments  to  be  in  error.  Our  calculations  show  that  the  self-energy 
correction  makes  the  0-2p  and  Ti-3d  bands  move  away  from  one  another  with  an 
accompanying  reduction  of  s,d  admixing,  in  agreement  with  the  latest  experiments. 
This  agreement  is  gratifying  since  this  method  makes  the  important  effect  of 
d-band  intrasite  correlations  and  the  formation  of  Hubbard  bands  accessible  to 
calculations . 
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The  near  normal  incidence  reflectivity  of  titanium  has  been  measured  for 
photon  energies  of  2  eV  to  25  eV  under  closely  controlled  surface  conditions. 
Measurements,  using  synchrotron  radiation  above  5  eV,  were  made  in  ultra  high 
vacuum  on  both  clean  and  in-situ  oxidized  surfaces  monitored  with  Auger 
spectroscopy.  The  reflectivity  data  were  Kramers-Kronig  analyzed  to  generate 
the  optical  constants.  Structure  in  Ej  appearing  in  the  oxidized  surface  but 
absent  in  the  clean  surface  is  related  to  the  oxygen  2p  and  2s  bands.  Some 
differences  in  the  oxide  structure  are  observed  between  oxide  films  grown  on 
titanium  exposure  to  oxygen  compared  to  exposure  to  an  oxygen  plasma. 
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INTRODUCTION 

The  optical  properties  of  titanium  and  its  oxides  have  not  been 

extensively  reported  in  the  literature  even  though  titanium  is  a  widely  used 

industrial  material  with  many  physical  and  chemical  properties  of  interest  inr 

materials  applications.  The  oxide  phases  have  a  variety  of  structures  and 

electrical  properties  which  range  from  metallic  behavior  in  TiO  to  insulating 

behavior  in  TiC^.  Suggestions  have  been  made  that  these  phase  changes 

accompany  oxide  growth  and  that  the  onset  of  passivation  in  the  oxide  films 

may  even  be  related  to  changes  in  phase. *  Most  of  the  optical  studies 
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attempted  have  been  on  metallic  titanium,  ’  although  the  present  work 
indicates  that  a  residual  oxide  was  present  in  these  previous  studies.  A 
systematic  study  of  oxide  growth  on  titanium  has  not  been  reported  in  which 
the  growth  of  the  oxide  was  monitored  with  optical  measurements. 

In  this  paper  we  report  measurements  of  the  optical  properties  of  oxides 
of  various  thicknesses  grown  on  titanium  under  a  variety  of  conditions.  The 
optical  properties  were  derived  from  reflectance  data  taken  in  the  range  2  to 
25  eV  with  some  data  extended  to  AO  eV.  The  data  were  Kramers-Kronig  analyzed 
and  the  values  for  and  z^  obtained.  Auger  spectra  were  taken  on  the  same 
surfaces  as  the  reflectance  measurements,  allowing  us  to  monitor  the  surface 
cleanliness.  The  analysis  of  the  Auger  lineshapes  was  critical  in  indicating 
the  possible  oxide  structures  formed.  An  analysis  based  on  the  reflectivity 
data  alone  would  have  been  insufficient  to  identify  the  oxide  structures. 

In  this  paper  we  compare  the  results  of  our  reflectivity  measurements 
with  the  results  of  other  measurements,  and  note  the  differences  among  the 
spectra.  Identification  of  the  structure  in  the  oxide  optical  constants  with 
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transitions  between  the  oxide  valence  and  conduction  bands  is  made  by 
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comparing  with  band  structure  calculations.  Analysis  of  our  Auger  spectra 
also  indicates  that  the  final  oxide  structure  for  our  oxygen  exposed 
polycrystalline  Ti  surface  but  not  our  oxygen  bombarded  surface  is  similar  to 
the  final  structure  observed  by  Shih  and  Jona^  for  an  oxygen  exposed  single 
crystal  surface. 


EXPERIMENTAL 


Near  normal  incidence  optical  reflectivity  and  Auger  spectroscopy 
measurements  were  performed  on  polycrystalline  titanium  samples.  All  of  the 
measurements  were  made  in  an  ultra  high  vacuum  system  with  a  base  pressure  of 
<1x10  ^  Torr.  Radiation  from  the  University  of  Wisconsin  Synchrotron 
Radiation  Center  was  used  as  the  light  source  for  photon  energies  greater  than 
5  eV  with  conventional  tungsten  and  discharge  lamps  used  for  the  low  energy 
end  of  the  spectrum.  The  sample  chamber,  illustrated  in  Figure  1,  is  designed 
such  that  optical  reflectivity  and  Auger  measurements  can  be  made  without 
repositioning  of  the  sample,  allowing  near  simultaneous  measurements.  A 
sputter  ion  gun  and  a  sample  heater  are  included  in  the  chamber  for  surface 
cleaning  and  oxidation  studies.  A  quadrupole  residual  gas  analyser  and  high 
purity  gas  manifold  are  also  attached  to  the  system.  Figure  2  shows  sche¬ 
matically  the  instrumentation  configuration.  Both  optical  and  Auger  measure¬ 
ments  are  performed  automatically  by  a  computer  based  data  acquisition  system 
described  elsewhere.^ 

Samples  12.7  mm  diameter  and  2  mm  thick  were  cut  from  commercially 
prepared  electron  beam  zone  .jlined  titanium  rod  having  a  nominal  purity  of 
99.972.  Major  impurities  included  oxygen  (,vl50  ppm)  and  Carbon  Ov  50  ppm). 

The  samples  were  mechanically  polished  with  a  final  polish  of  .05  micron 
alumina.  Following  the  mechanical  polish,  the  samples  were  placed  in  a 
trichloroethylene  vapor  degreaser  for  five  minutes,  washed  in  acetone,  and 
etched  for  approximately  5  seconds  in  a  solution  of  HF,  and  l^O. 

Following  the  etch  the  sample  was  washed  in  distilled  water  and  blown  dry. 

The  sample  was  immediately  mounted  and  placed  in  the  vacuum  chamber, 


pumped  down,  and  the  whole  system  baked  at  200°C  for  ^24  hours.  At  this  point 
the  Auger  spectra  show  strong  carbon  and  oxygen  lines,  while  the  titanium  line 
is  distorted  by  the  superposition  of  a  nitrogen  Auger  line.^  Heating  the 
sample  to  700°C  for  24  hours  removed  these  impurities;  however  large  Auger 

peaks  of  sulfur  and  chlorine  appeared,  presumably  segregating  on  the  surface r 
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from  the  bulk.  A  short  period  of  Ar  sputtering  (30  seconds,  lOpA/cm  at  3 

keV)  while  the  sample  was  being  heated  removed  these  contaminants,  leaving  the 

surface  clean  to  the  limits  of  detection  (approximately  1%  atomic).  LEED 

studies  performed  on  titanium  after  similar  cleaning  procedures  showed  an 
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ordered  surface.  The  optical  and  Auger  spectra,  after  sample  cooling,  would 
remain  stable  for  *vl  hour  in  the  MxlO  ^  Torr  vacuum.  A  short  heating  period 
of  30  minutes  to  600°C  followed  by  a  30  second  Ar  sputtering  at  600°C  was 
sufficient  to  reclean  the  sample.  All  reflectivity  measurements  on  the  clean 
sample  were  made  during  the  1  hour  "clean  period". 

A  series  of  optical  and  Auger  measurements  after  exposure  to  controlled 
amounts  of  oxygen  were  made.  The  oxygen  exposure  varied  from  5  Langmuir  (1 
Langmuir  =  10  ^  Torr  -  sec)  to  10,000  L,  all  at  room  temperature.  Subseq¬ 
uently,  measurements  were  made  on  the  sample  after  a  series  of  oxygen  ion 
bombardments  (energies  of  1  keV  and  2  keV).  Finally,  Ar  sputtering  and  sample 
heating  were  used  to  reclean  the  sample  in  order  to  check  the  reproducibility 
of  the  clean  surface  data. 


EXPERIMENTAL  RESULTS 

In  the  following  measurements  the  initial  surface  spectrum  is  that  of  the 
sample  when  first  placed  in  the  sample  chamber.  El  1 ipsometric  measurements  oi 
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Ti  under  similar  conditions  indicate  a  TiOj  surface  layer  of  about  85  A 

g 

thickness.  The  clean  surface  spectra  were  measured  following  the  cleaning 
procedure  previously  described.  The  oxidized  surface  spectra  occurred  after 
bombarding  the  sample  with  1  keV  oxygen  ions  at  15yA/cm  for  7.5  minutes. 

In  Figure  3  we  compare  the  reflectance  spectra  for  the  titanium  surfaces. 

We  attribute  the  structures  below  7  eV  to  transitions  within  the  d-band 

region  although  there  is  undoubtedly  some  s  and  p  mixing  in  these  bands.  This 

interpretation  is  consistent  with  that  of  Lynch  e^t  a_l.  ^  and  also  with  Ti 

I'll  III  x-ray  emission  and  absorption  data^  which  give  an  indication  of  the 

shapes  and  relative  energy  positions  of  the  valence  and  conduction  bands, 

respectively.  We  attribute  the  peak  around  10  eV  and  the  structure  beyond  20 

eV  in  the  oxidized  and  initial  surface  spectrum  to  transition  from  the  oxygen 

2p  and  2s  bands,  respectively,  although  the  width  of  the  latter  structure 

indicates  transition  from  the  valence  to  higher  bands  may  also  be  present. 

These  peak  assignments  are  consistent  with  energy  band  calculations  which 

indicate  transitions  involving  the  oxygen  2s  and  2p  bands  should  occur  at 

these  energies.  In  addition,  the  separation  of  16  eV  between  the  two  peaks 

in  the  initial  surface  spectrum  is  just  the  free-atom  value  and  has  been  seen 
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in  the  XPS  spectrum  of  TiC^.  It  is  to  be  noted  that  these  two  peaks  are  at 
lower  energy  in  the  initial  surface  spectrum  as  compared  to  the  oxide  spect¬ 
rum.  We  believe  the  spectra  differ  because  a  TiO-like  phase  exists  on  the 
oxidized  surface  whereas  reduced  TiC>2  exists  on  the  initial  surface.  These 
oxides  are  not  stoichiometric,  and  the  air-exposed  surface  may  even  be 
hydrated.  However,  our  viewpoint  is  corroborated  by  the  Auger  data  which  show 
an  oxygen  peak-to-peak  height  twice  as  large  for  the  initial  surface  as  for 
the  oxidized  surface,  and  also  show  lineshapes  indicative  of  these  oxide 
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phases.  In  addition,  there  is  a  reduction  in  the  height  of  the  d-band 

structure  between  5-7  eV  relative  to  the  oxygen  2p  band  in  the  initial  surface 

Reflectance  as  compared  to  the  oxidized  surface.  Reduced  Ti(>2  should  have  a 

qujch  lower  occupation  in  the  d-band  region  than  TiO.  Recent  UPS  measurements 
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of  TiO  exposed  to  oxygen  show  that  there  is  indeed  a  dramatic  reduction  of  * 
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the  d-band  structure  under  oxygen  exposure.  In  these  results  the  Fermi 
level  also  moves  about  1  eV  closer  to  the  valence  band  as  the  sample  is 
oxidized  with  a  position  at  maximum  oxygen  exposure  the  same  as  for  stoich¬ 
iometric  TiO^.  This  movement  is  consistent  with  the  change  in  position  of  the 
oxygen  2p  band  displayed  in  Figure  4. 

For  comparison  we  present  the  spectrum  of  Lynch  et  al.^  and  our  cleanest 
surface  spectrum  in  Figure  4.  The  peaks  around  10  eV  and  25  eV  appearing  in 
Lynch's  data  disappear  in  ours,  although  a  shoulder  at  about  13  eV  remains. 

The  Lynch  spectrum  falls  somewhere  between  our  initial  surface  and  oxidized 
surface  spectra. 

Figure  5  shows  the  reflectance  spectra  after  successive  oxygen  bombard¬ 
ments  with  the  build-up  of  the  oxygen  peak  at  about  10  eV,  the  movement  to 
lower  energy  and  build-up  of  the  oxygen  peak  above  20  eV,  and  the  reduction  of 
the  d-band  structure  between  5-7  eV.  Note  that  although  the  bombardment 
builds  an  oxide  much  thicker  than  that  on  the  initial  surface,  none  of  the 
spectra  match  the  initial  surface  spectrum.  In  particular  the  peak  at  10  eV 
remains  fixed  in  energy  even  for  the  thick  oxides,  and  does  not  move  to  lower 
energy  as  for  the  initial  surface.  However,  the  peak  above  20  eV  does  not 
move  to  lower  energy  and  the  d-band  structure  is  reduced,  these  features  being 
comparable  to  the  initial  surface  structure.  It  is  possible  that  during  the 
oxygen  bombardment  regions  of  different  oxide  structure  are  formed  on  the 
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surface.  The  predominant  structure  might  be  TiO-like,  but  with  regions  having 
TiO^  structure.  The  reflectance  spectrum  would  then  show  a  combination  of  the 
features  of  these  oxides.  The  fact  that  a  sub-stoichiometric  oxide  may  con¬ 
tinue  to  grow  under  oxygen  bombardment  to  a  thickness  much  greater  than  the 
initial  surface  thickness  indicates  different  conditions  of  exposure  can 
produce  quite  different  oxide  structures.  This  result  may  be  of  consequence 
in  the  study  of  the  mechanisms  of  oxide  growth  and  in  the  search  for  means  to 
control  oxide  growth  rates. 

We  must  analyze  the  Auger  lineshapes  in  detail  to  obtain  further  infor¬ 
mation  about  the  surface  oxide  structure.  Our  Auger  spectrometer  has  an 
energy  resolution  AE/E  of  about  .6%  so  that  these  measurements  in  the  range 
from  300  to  500  eV  have  a  resolution  of  2  to  3  eV.  The  Auger  spectra  in  Fig. 

6  are  for  the  surface  depicted  in  the  reflectance  spectra  of  Fig.  3.  The 

figure  shows  the  Ti  ^  and  0  KLL  structure.  The  clean  surface  Ti  L^  jjj 
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spectrum  is  identical  to  those  of  Solomon  and  Baun  for  a  freshly  evaporated 
film  and  Shih  and  Jona^  for  the  Ti(001)  surface.  The  initial  surface  and 
oxidized  surface  spectra  have  peak  to  peak  ratios  between  the  dominant  Ti  LMM 
structure  and  the  0  KLL  structure  which  when  adjusted  for  transition  probabil¬ 
ities  are  2:1  and  1:1,  respectively.  This  result  suggests  a  TiOj  structure  on 

the  initial  surface  and  a  TiO  structure  on  the  oxidized  surface.  The  initial 
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surface  spectrum  is  nearly  the  same  as  the  TiOj  Auger  spectrum  shown  in  Fig. 
4  of  Ref.  14  except  that  the  shoulder  at  about  420  eV  is  more  pronounced  in 
our  spectrum.  The  most  prominent  difference  between  our  oxidized  surface 
spectrum  and  the  TiO  spectrum  in  Fig.  4  of  Ref.  14  is  the  absence  in  our 
spectrum  of  the  shoulder  feature  at  about  410  eV  in  the  Ti  LMM  structure. 


We  studied  the  effects  of  oxidation  under  different  conditions  on  the 
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Auger  spectra.  The  above  oxygen  bombarded  surface  corresponds  to  about  500  L 
oxygen  exposure.  At  about  100  L  exposure  to  oxygen  gas  at  room  temperature, 
the  Auger  spectrum  of  Shih  and  Jona^  for  Ti(001)  was  already  identical  to  the 
TiO  spectrum.  Further  oxygen  bombardment  for  several  more  minutes  produced 
the  changes  in  the  reflectivity  noted  in  Fig.  5,  but  the  Auger  spectrum  did  * 
not  change  much  and  in  particular  the  Ti/0  peak  ratio  remained  the  same. 

We  also  studied  the  effects  of  oxidation  at  room  temperature  cn  our 
polycrystalline  Ti  sample  .  In  this  case  the  pattern  of  change  in  the  Ti 
L , ,  spectrum  observed  by  Shih  and  Jona^  was  also  observed  by  us,  but  the 

changes  were  much  slower  with  respect  to  oxygen  exposure.  Shown  in  Fig.  7  are 
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the  spectra  at  500  L  and  10  L  oxygen  exposure.  At  500  L  one  can  clearly  see 
the  feature  at  A10  eV,  but  it  is  still  less  pronounced  than  in  the  Shih  and 

4 

Jona  spectrum  at  100  L.  After  10  L  exposure,  our  spectrum  is  identical  with 
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that  of  Shih  and  Jona  and  also  with  the  TiO  spectrum  .  At  500  L  the 
reflectance  spectrum  is  nearly  the  same  as  our  oxidized  surface  spectrum  in 
Fig.  3;  and,  unlike  the  oxygen  bombarded  case,  the  reflectance  remains  about 
the  same  for  greater  oxygen  exposures. 

KRAMERS  KR0NIG  ANALYSIS 

A  Kramers-Kronig  analysis  was  used  to  derive  the  phase  6  of  the  complex 
reflectance  re1**  from  the  reflectivity  data  and  thus  to  derive  the  complex 
refractive  index  and  the  complex  dielectric  function.  For  energies  below  2 
eV,  the  same  reflectance  data  segment,  taken  from  measurements  by  Lynch 
et  al.,***  was  used  throughout.  As  we  discuss  below,  this  low  energy  approx¬ 


imation  did  not  affect  the  location  of  structure  in  optical  constants  derived 
from  the  KK  analysis  in  the  region  of  2-25  eV. 
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Since  an  infinite  integral  of  the  reflectivity  over  frequency  is 
necessary  to  obtain  the  phase  information,  a  high  energy  extrapolation  of  the 
data  was  needed.  We  used  a  standard  power  law  extrapolation  of  the  form 
R(E)  ■  RqE  for  the  reflectivity.  Such  extrapolations  are  usually  not 
accurate  because  of  high  energy  structure  present  in  most  materials,  unless  r 
the  exponent  6  can  be  fitted  to  give  a  known  quantity,  such  as  an  indepen¬ 
dently  obtained  optical  constant,  at  a  given  energy.^  However,  we  have  been 
able  to  obtain  consistent  results  using  a  method  which  requires  no  indepen¬ 
dently  obtained  optical  parameters.  We  fit  the  parameter  3  such  that  one  of 


the  sum  rules 
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is  satisfied.  The  latter  constraint,  Eq.  (2),  depends  less  on  high  frequency 
structure  than  does  the  former^,  and  thus  we  use  Eq.  (2)  to  fit  the 
extrapolated  curve.  For  the  data  considered  here,  Eq.  (1)  also  has  a 
reasonably  small  value  in  all  cases. 

The  various  segments  of  the  reflectivity  spectrum  from  the  discharge 
source  and  from  the  synchrotron  had  to  be  scaled  to  provide  a  smooth  curve. 

In  order  to  test  the  effect  on  the  optical  constants  of  this  scaling  and  of 
the  extrapolation  procedure  used,  we  calculated  the  optical  constants  of  the 
oxidized  titanium  surface  for  two  extreme  cases  of  scaling.  In  both  cases  the 
parameter  8  was  fit  to  satisfy  Eq.  (2).  Figure  8  compares  the  joint  density 
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of  states  for  these  two  scalings.  It  is  seen  that  the  spectral  shapes  compare 
well  up  to  about  18  eV  with  peak  shifts  of  no  more  than  a  few  percent,  but 
that  the  final  peak  at  22-25  eV  shifts  by  3-4  eV,  as  is  to  be  expected  since 
this  peak  is  at  the  edge  of  the  data.  Since  the  shapes  of  our  reflectivity 
spectra  are  given  much  more  accurately  than  their  magnitudes,  the  positions  ' 
and  relative  magnitudes  of  the  peaks  in  the  optical  constants  are  reliable, 
with  some  shifting  possible  for  the  structure  above  20  eV.  Note  tl.jt  in 
Figure  8  the  scaling  has  not  affected  the  spectral  shape  below  2  eV. 


ANALYSIS  OF  OPTICAL  CONSTANTS 

The  complex  refractive  index  n-ik,  derived  from  our  reflectances  and 
phase  results,  can  be  used  to  construct  the  complex  dielectric  function  via 
the  expression 

e1  -  ie2  -  (n-ik)2 

The  quantity  z^  gives  the  absorption  spectrum  of  the  system  and  therefore 
describes  the  excitation  structure.  If  one  neglects  final  state  effects 
introduced  by  the  interaction  of  the  excited  electron  with  its  hole  and  also 
neglects  any  intraband  contributions  so  that  only  direct,  interband 
transitions  contribute,  then  EjCw)  has  c^e  f°rm 

e7(u>)  %  4  1  E  „£)“«>  K  a£>|2  d3k  <3> 

2  n,s  BZ  n,S  n,S 

The  subscripts  n  and  s  refer  to  filled  and  unfilled  bands  respectively,  and 

u>  (k)  (k)  -  0)  (t).  The  momentum  matrix  element  M  (k)  between  states 

n,s  n  s  n,s 

in  bands  n  and  s  is  nearly  constant  for  bands  with  little  dispersion.  Then 

e-(u>)  N.(u>)  (A) 

*  0)*  J 

where  Nj(w)  is  the  joint  density  of  states. 
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The  KK  analyzed  spectra  of  and  for  the  clean,  the  oxidized,  and  the 
initial  titanium  surfaces  are  shown  in  Figure  9.  The  changes  in  structure 
among  the  £2  spectra  are  in  most  cases  similar  to  the  changes  in  the  reflect** 
ivity  spectra  in  Figure  4.  Note,  however,  that  the  peak  position  at  about  10 
eV  in  the  oxidized  spectrum  is  closer  to  the  position  in  the  initial  1 

spectrum  than  are  the  peak  positions  in  the  corresponding  reflectivity 
spectra.  The  struct^  is  also  less  pronounced  in  the  £2  spectra  than  the 
corresponding  structure  in  the  reflectivity  data. 

To  gain  a  more  detailed  understanding  of  the  oxide  structure,  we  may 

compare  our  results  with  the  densities  of  states  calculated  by  Schoen  and 

Denker  (SD).^  These  calculations  take  the  defect  character  of  the  oxide  into 

account,  allowing  for  a  15%  vacancy  concentration  and  the  non-stoichiometry  of 

the  oxide,  and  the  valence  bands  are  in  substantial  agreement  with  the  most 
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recent  XPS  measurements  on  Ti01+x.  Though  we  do  not  have  the  joint  den¬ 
sities  of  states  from  these  calculations,  we  may  note  from  Eqs.  (3)  and  (4) 
above  that  structures  in  the  density  of  states  which  are  at  different  energies 
but  are  from  the  same  regions  in  the  Brillouin  zone  can  be  related  to 
structure  in  the  optical  spectrum.  Upon  studying  Fig.  1  in  SD,  showing  the 
band  structure  and  density  of  states  for  stoichiometric  TiO  with  15% 
vacancies ,  we  see  that  indeed  the  prominent  structures  there  can  be  linked  to 
the  structure  in  the  joint  density  of  states  for  the  oxidized  surface,  the 
topmost  curve  in  Fig.  8.  The  0-2p  density  of  states  will  link  strongly  to  the 
mostly  d-like  first  peak  in  the  conduction  band,  the  peak  regions  in  the  two 
densities  coming  mostly  from  parallel  bands  in  the  Brillouin  zone  with  the 
peaks  separated  by  about  9  eV,  in  agreement  with  our  assignment  of  the 
structure  in  Fig.  8.  The  0-2s  band  will  have  a  weaker  transition  to  the  s-p 
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hybridized  bands  centered  about  1.25Ryd  in  the  density  of  states  giving  a 
separation  of  about  26  eV,  again  in  agreement  with  experiment,  although  the 
breadth  of  the  experimental  structure  may  indicate  transitions  from  the 
valence  bands  to  higher  conduction  bands  also  contribute.  The  characteristic 
d~band  double  peak  structure  at  about  5-6  eV  in  Fig.  8  may  be  due  to  tran-  * 
sitions  to  the  double-peaked  hybridized  s-p  density  of  states  in  Fig.  7  of  SO, 
or  it  may  be  due  to  structure  in  the  valence  d-band  density  of  states  itself, 
as  indicated  by  the  calculation  of  Em  and  Swi tendick. * ^  The  results  of  SD 
for  other  compositions  of  TiOj+x  also  support  our  assignments  of  transitions, 
although  the  transition  energies  would  vary  somewhat. 

Comparison  with  the  energy  band  calculations  for  TiO  allows  assignments 
of  the  peaks  in  the  optical  spectrum,  but  the  calculations  in  themselves  do 
not  enable  us  to  identify  the  oxide  structure  as  TiO.  Comparison  with  other 
band  structures  for  different  oxide  compositions  would  probably  allow  similar 
assignments  for  the  optical  constants.  However,  the  identification  of  the 
peaks  in  the  5-6  eV  region  of  the  oxidized  spectrum  as  d-transitions  and 
there  subsequent  decrease  under  further  oxygen  bombardment  or  air  exposure 
indicates  much  further  oxidation  occurs  in  these  cases.  Still,  the  Auger 
lineshapes  are  necessary  to  make  clear  identifications  of  the  oxide  strucure. 

CONCLUSIONS 

The  reflectivities  of  titanium  surfaces  undergo  significant  changes  as 
oxides  are  formed.  Careful  cleaning  monitored  by  Auger  spectroscopy  is 
necessary  to  insure  a  clean  Ti  surface,  and  several  previous  Ti  reflectivity 
measurements  show  significant  structure  due  to  residual  oxides.  Cycles  of 
oxygen  bombardment  followed  by  Argon  sputtering  yield  reproducible 
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reflectivity  spectra.  The  previously  unidentified  peaks  at  about  8-10  eV  and 
26  eV  in  the  reflectivity  spectra  are  most  likely  due  to  transition  from 
oxygen  2p  and  2s  bands,  respectively.  Band  structure  calculations,  x-ray 
emission  and  absorption  data,  and  x-ray  photoelectron  data  which  show  similar 
band  separations  for  several  other  transition  metal  oxides  support  this 
conclusion. 

The  oxide  produced  initially  by  oxygen  bombardment  on  the  clean  surface 
seems  to  be  TiO-like  with  further  bombardment  producing  a  trend  toward  greater 
oxidation  or  perhaps  a  mixed  oxide  state,  but  the  air  exposed  initial  surface 
seems  to  be  TiO^.  Oxygen  exposure  of  the  clean  surface  at  room  temperature 
eventually  produces  a  stable  TiO  structure,  but  the  structure  evolves  at  much 
greater  oxygen  exposures  than  for  the  single  crystal.^ 

The  changes  in  the  reflectivity  data  which  accompany  these  different 
conditions  of  oxygen  exposure  are  consistent  with  the  assumed  oxide 
structures,  but  the  reflectivity  in  itself  is  not  sufficient  to  identify  the 
oxide  structures.  Comparisons  of  Auger  peak  heights  and  especially  invest¬ 
igations  of  Auger  lineshapes  are  necessary  to  identify  oxide  structures.  On 
the  other  hand,  certain  trends  in  the  reflectivity  data,  such  as  shown  in  Fig. 
5,  are  not  readily  apparent  in  the  Auger  data.  Both  types  of  measurement  were 
necessary  in  this  work  to  identify  important  trends  in  the  oxide  growth. 

Other  investigators  have  seen  changes  in  oxide  structure  under  different 

1 8 

conditions  of  exposure.  Recent  measurements  by  Garner  et  al. ,  of  the 
initial  oxidation  stages  of  the  Si(lll)  surface  indicate  that  different  oxide 
structures  form  depending  on  whether  the  oxygen  gas  is  excited  or  not.  The 
SiOj  structure  occurs  only  for  exposure  to  excited  oxygen.  Undoubtedly  the 
condition  of  oxygen  bombardment  in  our  measurements  prevents  the  surface  from 

.  \ 
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stabilizing,  and  our  air  exposed  oxide  surface  may  have  a  more  complicated 
structure,  such  as  a  hydrated  structure  which  would  be  difficult  to  detect. 
Certainly  the  differences  among  these  oxides  may  be  quite  important  in  the 
mechanisms  of  oxide  growth  under  various  conditions,  and  we  are  planning 
further  investigations  along  these  lines. 
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Figure  Captions 

Cutaway  view  of  ultrahigh  vacuum  chamber  with  apparatus  for 
Auger  and  optical  reflectivity  measurements. 

Instrumentation  configuration  for  the  experimental  chamber. 

The  ellipsometer  instrumentation  shown  is  under  development 
and  was  not  used  in  this  work. 

Reflectance  spectra  for  various  Ti  surfaces. 

Comparison  of  our  cleanest  surface  spectrum  with  the  measurements 
of  Lynch,  Olson,  and  Weaver.10 

Reflectivity  of  Ti  under  oxygen  bombardment.  The  top  curve  is  the 
cleanest  surface,  and  the  curves  below  are  after  successive  oxygen 
bombardments.  Each  set  of  data  is  offset  for  clarity  with  the 
bracket  indicating  a  5  percent  change  in  reflectivity. 

Auger  spectra  taken  just  prior  to  the  optical  measurements.  The 
positive  excursion  of  the  first  major  Ti  core-core  Auger  line  in 
the  initial  and  oxidized  surface  spectra  is  due  to  a  coincident 
nitrogen  Auger  line.  Note  that  the  ratio  of  the  0  to  Ti  peak- 
to-peak  heights  is  twice  as  large  in  the  initial  surface  spectrum 
as  in  the  oxidized  surface  spectrum. 

4 

Ti  LMM  Auger  spectra  at  500  L  and  10  L  oxygen  exposure.  The  latter 
spectrum  is  identical  with  that  for  TiO.14 

Joint  density  of  states  for  the  Kramers-Kronig  analyzed  reflectivity 
data  using  two  extreme  cases  of  scaling  of  the  raw  reflectivity  data 

Optical  constants  obtained  from  the  Kramers-Kronig  analysis  of  the 
reflectivity  data.  Each  curve  is  offset  as  indicated.  (a)  The 
real  part  of  the  dielectric  constant.  (b)  The  imaginary  part  of 
the  dielectric  constant. 
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